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Abstract 

In 1960 Schwinger [J. Schwinger, Proc.Natl.Acad.Sci. 46 (1960) 570- 
579] proposed the algorithm for factorization of unitary operators in 
the finite M dimensional Hilbert space according to a coprime decom- 
position of M. Using a special permutation operator A we generalize 
the Schwinger factorization to every decomposition of M. We obtain 
the factorized pairs of unitary operators and show that they obey the 
same commutation relations as Schwinger's. We apply the new factor- 
ization to two problems. First, we show how to generate two A;g-like 
mutually unbiased bases for any composite dimension. Then, using a 
Harper- like Hamiltonian model in the finite dimension M = M1M2, 
we show how to design a physical system with M\ energy levels, each 
having degeneracy M%. 



1 Introduction 

A finite phase space of dimension M, where coordinate and momentum have 
M possible values, is a frequent component of various physical and mathe- 
matical problems. Fast Fourier Transform (FFT) [TJI2], Schwinger factoriza- 
tion of unitary operators p], generation of kq bases and finite dimensional 
Harper-like Hamiltonians [I] are the problems related to finite phase space 
that will be considered in this paper. A recent review of various quantum 
systems with finite Hilbert space can be found in ref. [5]. 



Originally studied by Weyl [5], the finite dimensional Hilbert space was 
systematized by Schwinger in terms of "Unitary Operator Bases" |3J. Schwinger 
considered a M-dimensional physical system. Such a Hilbert space can be 
achieved by application of the following boundary conditions on the wave 
function ip{x) and its Fourier transform ^(p) (ref. [7]): 

27r/i 

i/j(x) = if)(x + Mc), <5/(p) = ^(p^ ), (1.1) 

c 

where c is a length unit. In what follows, we will assume c = 1. As a 
consequence of the above boundary conditions x and p have a finite discrete 
spectrum of eigenvalues: 

x = 0,l,...,M-l; p =^.{0,l,...,M-l}. (1.2) 
Using unitary operators U and V (ref. [1] with c = 1): 



U = e i£ %, 
V = e&, 



1.3) 



the complete orthogonal operator basis of M 2 operators can be defined as 
®- 

U k V n ; k,n = 0,1,..., M- 1. (1.4) 
The above operators have the commutation relation: 

V n U k = U k V n e^ nk . (1.5) 

For a coprime decomposition of M = MiM 2 , using the Fermat-Euler theo- 
rem, Schwinger showed how to factorize the unitary operators. The Fermat- 
Euler theorem states that if Mi and M 2 are coprime, then there exist unique 
iVi and A^2 such that: 

M ± N 2 = 1 (mod M 2 ), M 2 N 1 = 1 (mod Mi). (1.6) 
Therefore, the two pairs of unitary operators defined as: 

Ui = U M \ \ U 2 = U M \ 



v 1 = v h,hN \ \v 2 = v A ' hN \ ( LT ) 

behave as independent complementary operators of factorized dimensions Mi 
and M 2 . The respective commutation relations are: 

y7«U*i = U-*V^e^ n * k \ i = 1, 2; (1.8) 
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\r i''; j c/\r ■< / j- i,j = 1,2. 



(1.9) 



Each unitary operator on the M dimensional space (Eq. I1.4p can be 
considered as two operators from the factorized dimensions Mi and M 2 . This 
is due to the one-to-one "Sino-Ruritanian" correspondences [8J: 



Therefore, for every power n of the operator V we can find the unique rep- 
resentation by the factorized unitary operators V\ and V 2 . The appropri- 
ate powers rii and n 2 of the factorized operators V\ and V 2 are determined 
by the first "Sino-Ruritanian" correspondence (Eq. (11.101) ). Similarly, the 
correspondence between U and (1/1,1/2) is determined by the second "Sino- 
Ruritanian" correspondence (Eq. (ll.lOp ). Another recent factorization con- 
struction based on the Chinese Remainder Theorem (CRT) can be found in 



After we have obtained factorization of the M dimensional Hilbert space 
into its coprime sub-dimensions Mi and M 2 , we can apply it to the kq bases 
generation and the Harper-like Hamiltonian model. Let us first consider the 
kq bases generation. The factorized operators from Eq. ( 11. 7p can be used 
for generation of the following two pairs of operators (note that V 2 Ml = V Ml 



where the dimension M = MiM 2 and a = M 1; b = M 2 (according to the 
notation of ref. jl] with c=l). Hence, by employing all possible powers, 
each pair of operators (a) and (b) forms a complete set of M commuting 
operators and thus generates an alternative kq basis for treatment of the 
M-dimensional Hilbert space. We have two such bases: 



(«) I*, ?> = 7fe E^o" 1 + sa), (b) \K, Q) = -fa e iKth \Q + tb) 



n = n l M 2 N l + n 2 M 1 N 2 (modM) 
k = k 1 M 2 + k 2 M 1 (modM). 



(1.10) 



ref. 0. 



and V X M2 = V M2 ): 




(1.11) 




(1.12) 
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The unique property of the kq bases is that they are eigenf unctions of both 
space and momentum displacement operators. These functions have partial 
knowledge about both position and momentum, whose precise simultaneous 
knowledge is limited by the non-commutation of the corresponding operators. 
In the case of dimension M = M\M 2 factorizable to coprime numbers M\ 
and M 2 , the two kq bases (a) and (b) are Mutually Unbiased Bases (MUB) 
|3]. The MUB property of the bases means that if the physical system is 
found in one of the states of one MUB (for example set (a)), then it has 
equal probabilities to be in all the states of the other MUB (set (b) in our 
example). Mathematically, the mutual unbiasedness of the two kq bases 
means the following equality: \(k, q\K, Q)\ 2 = jj. For non-coprime M\ and 
M 2 the MUB property is violated. For example, if Mi = m\T and M 2 = m 2 r 
have a common multiple r, we have: 

(k, q\K,Q) = -= e- lksa e lKt \sm 1 r + q\tm 2 r + Q) . (1.13) 

t,s 

The product (srriir + q\tm,2r + Q) equals unity for the solution of the following 
modular equation: 

srriir + q — tm 2 r — Q = (modM). (1-14) 

Following ref. [H] (p. 45 theorem 'd'), the above equation can be taken 
modulo r: 



q = Q(modr). (1-15) 

Therefore, for q = and Qq = 1 (which is always possible according to the 
ranges of values Eqs. (I1.12|) ) we have \(k, qo\K, Qo)\ 2 = 7^ ^. 

To complete the introduction to kq MUB we note their quasi-periodic 
properties: 



{a)\k+^,q) = \k,q), \k,q + M 1 ) = e~ ika \k, q), 
{b)\K + £,Q) = \K,Q), \K,Q + M 2 )=e-* Kb \K,Q). 

Now, let us consider Harper-like Hamiltonians. They are defined as 
Hamiltonians of one degree of freedom periodic both in coordinate and mo- 
mentum [10]. For our discussion we are interested in the use of Harper-like 
Hamiltonians for the energy spectra design considered in ref. The energy 
spectra design is a direct consequence of the factorization of the M = M\M 2 
- dimensional Hilbert space to coprime constituents M\ and M 2 . In the origi- 
nal version (ref. |4J) one considered a Harper-like Hamiltonian H[T(b), r(— )] 
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which is a function of the two operators T(b) and r(— ) from Eq. fll.lip . It 
is important that the Hamiltonian H[T(b), r(— )] is a function of the oper- 
ators V\ and Ui (due to T(b) = V^ 2 and r(— ) = U\)] in such a case only 
the Mi - dimensional subspace is affected by the Hamiltonian. The M 2 - 
dimensional subspace is untouched by the Hamiltonian. Hence, considering 
H[T(b), we expect to obtain Mi energy levels (with a spectrum de- 

termined by the details of the Hamiltonian) each of which is degenerate M 2 
times. 

The aim of this paper is first to extend the Schwinger factorization to 
non-coprime M\ and M 2 . Then the other two applications, the A;g-like MUB 
generation and the energy spectra design by Harper-like Hamiltonian, are 
extended correspondingly. For that purpose, in section 2, we define the 
permutation operator A, based on the previous study by Cooley and Tukey 
of Fast Fourier Transform (FFT) [2]. Using the operator A we obtain pairs 
of unitary operators, which have commutation relations as in Eqs. (II .8} II .9p . 
In section 3 we use the new factorized unitary pairs to generate two kq-\ike 
MUB. New quasi-periodicity properties are obtained in one of the bases. In 
section 4 we apply the new factorization to the energy spectra design using 
Harper-like Hamiltonians without any restriction on the factors Mi and M 2 
of the dimension M = M\M 2 . Section 5 includes a discussion and summary. 

2 Factorization of unitary operators using the 
permutation operator A 

To define the permutation operator A we start by recalling the Division 
Algorithm Theorem (DAT) from number theory. 

The theorem states (ref. [9] page 2 or ref. |TT] page 3) that for any integer 
numbers D and d with d > 0, there exists a unique pair of integer numbers 
q and r satisfying the following conditions: 



Consider the special case of positive integer D in the range [0,1,..., M\M 2 — 
1] and positive d = M 2 . In this case there is a unique pair of integers q and 
r satisfying the following conditions: 



(a) D = d ■ q + r, 
(6) < r < d. 



(2.1) 



(a) q G [0,1,..., 

(b) re [0,1,..., 



Mi-1] 
M 2 -l] 



(2.2) 



(c) D = d ■ q + r. 
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For our discussion, this DAT based special representation of the numbers 
modulo M = M\M 2 is the crucial component. 

In 1965 Cooley and Tukey [2] introduced an FFT algorithm not limited 
to coprime factorization of M — MiM 2 . They used two complementary DAT 
based representations for the x and p variable indices: 



n = n 1 M 2 + n 2 , . , 

k = k 1 + k 2 M 1 , 

which enabled them to simplify the Discrete Fourier Transform (DFT) cal- 
culation (wm = e"^): 



„ _ _J_ V^M-l nk _ _J_ v-M-1 (n 1 M2+n 2 )(k 1 +k2M 1 ) _ 

(2.4) 

1 n 2 (fci+fc 2 Mi) l rnfcjMa 

^ 2^112=0 M y/M{ Z>ni=0 x nUJ M 

The two summations in the last line of the above formula require M- Yli=i Mi 
operations instead of M 2 operations by direct calculation [2]. 

Following Cooley and Tukey, we define a permutation operator A, which 
acts in the finite M - dimensional Hilbert space: 

Afi-l M 2 -l 

A= \X2 + M 2 x 1 )(x l + M 1 x 2 \. (2.5) 

xi=0 X2=0 

For the construction of the operator A we used two DAT based representa- 
tions, as in Eq. (12. 3ft . applied to the coordinate states \x). In the coordinate 
representation our operator is equal to the stride permutation matrix widely 
used in signal processing |12| . For a simple illustration, let us consider the 
example of dimension M = 6, where Mi = 2 and M 2 = 3. The table of 
correspondence between the numbers x = 0, 1, 2, 3, 4, 5 and pairs of numbers 
(xi = 0, 1; x 2 = 0, 1, 2) according to the rule x = X\ + M\X 2 is: 



X 





1 


2 


3 


4 


5 


Xi 





1 





1 





1 


x 2 








1 


1 


2 


2 



Table 1: DAT representation of x = (xi,x 2 ) 
With the rule x = x x M 2 + x 2 we have another table : 
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X 





1 


2 


3 


4 


5 


1 

x x 











1 


1 


1 


x 2 





1 


2 





1 


2 



Table 2: DAT representation of 



Consequently, the permutation matrix A corresponding to Eq. 
ing the standard basis for \x)) is: 



A 



us- 



/ 1 














o\ 








1 























1 








1 























1 








\o 














1 / 



(2.6) 



Also, it can be written in a compact way as A = (0)(1, 3, 4, 2)(5). This means 
that A leaves the coordinate states |0) and |5) unchanged, |1) goes into |3), 
|3) goes into |4), |4) goes into |2) and |2) goes into |1). We note that the 
operator A is unitary: 



AA ] = I. (2.7) 

With the permutation operator A at hand, we can define the new fac- 
torization. To do this we use the pairs of operators from Eq. ( II. lip , whose 
definition captures a general factorization of M = M1M2. We modify the 
(a) set of operators by the permutation operator A of Eq. (I2.5p . and for 
convenience relabel all the operators: 

(n >) f Ui = r'(f ) = At(*)A1 = AU^A; f U 2 = r(f ) = e«¥ = U M ^ 
{a) \ V 2 = T'{a) = AT(a)A^ = AV Ml A^\ K ' \V X = T(b) = e^ b = V M \ 

(2.8) 

The tilde denotes the new version of operators. The (a') and (b) sets of 
operators replace the Schwinger operators of Eq. (II. 7p . As we will show 
shortly, they obey all the commutation relations Eqs. (11.81 II. 9p of factorized 
operators. Therefore, the (a') and (b) sets of operators define the new factor- 
ization, not restricted to coprime decomposition. The commutation relation 
Eq. (II. 9p is fulfilled by the tilde operators (Eq. (12. 8p ) due to the unitarity 
property of A. For the commutation relation Eq. (11.81) we first calculate the 
operation of AU M2 A^ and A^U Ml A on coordinate states: 
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AU M2 A*\x) = AU M2 A^\ Xl M 2 + x 2 ) = AU M2 \x x + M x x 2 ) 
= e^ Xl A\xi + M x x 2 ) = e^ Xl \x 1 M 2 + x 2 ); 



A^U Ml A\x) = A^U Ml A\x\ + M x x 2 ) = A^U Ml \x[M 2 + x' 2 ) 
= M- x% A^\x x M 2 + 4) = e^ X2 \x' x + Mi4), 



(2.9) 



(2.10) 



where the only difference in the above calculations is that we used different 
DAT based representations for the x values. Using the above expressions one 
can prove: 

V^U^lx) = V M2ni AU Ahkl A^\ Xl M 2 + x 2 ) = e^ klxl V M ^\xiM 2 + x 2 ) = 

= e^ klXl \(x 1 -n 1 )M 2 + x 2 ), 

(2.11) 

whereas applying U^V™ 1 we get: 

U^V^\x) = AU M * kl Aty M ™\ Xl M 2 + x 2 ) = AU Makl Ai\(xi - n x )M 2 + x 2 ) = 

= ( M k ^- ni) \{x l -n l )M 2 + x 2 ). 

(2.12) 

Similar results can be shown for the operators V 2 2 and U 2 2 . Summarizing 
the results, the commutation relation (Eq. (11.81) ) is fulfilled: 

V7*ffki = (/^"-et"-'; i = 1,2. (2.13) 

Therefore, using the permutation operator A we obtained the new factoriza- 
tion of the unitary operators, which is not limited to coprime decomposition 
of M. Here we used the permutation operator A for the transformation of 
the (a) set of operators to the new (a') set for the factorization. Obviously 
we could have applied the transformation to the (b) set, which would have 
also enabled the factorization. 

In the particular case of Mi = M 2 (a=b), the two kq bases Eq. (11. 12ft 
are identical, and so are the two (a) and (b) sets of operators in Eq. (II. lip , 
and the permutation operator A from Eq. (12. 5p satisfies A 2 = I. In this case 
we have only one set of kq - operators and only one \k, q) basis. Application 
of the permutation operator A to that set of operators defines the tilde set, 
which obeys the proper commutation relations with the original set (Eqs. 
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II. S\ ll.9p . Their respective eigenstates (obtained by applying A to the unique 
\k,q) basis) are MUB with respect to the original \k,q) states. (See also the 
next example and the treatment of Harper - like Hamiltonians for M = 4 = 2 2 
in section 4). 

The permutation operator A, based on the analogy to the Cooley and 
Tukey FFT, solves the unitary operator factorization. To acquire some phys- 
ical intuition about the operator A, let us consider the example of dimension 
M = 4, where M\ = Mi = 2. In this case, the operators from Eq. fll.lip in 
the coordinate representation may be presented as: 



(«) 



U 2 



( 1 








^ 









1 


M 





-1 








, v 2 = 











1 








1 





1 











\o 








-1 J 




Vo 


1 





o / 



(2.14) 



The sets (a) and (b) of operators in Eq. fll.lip are identical in our example. 
To get the second set of factorized operators we write the (a') set from Eq. 



(a) I 



AU 2 A^ 



( 1 








^ 





1 














-1 





\o 








-1 > 



AV 2 A^ 





i 





°\ 


1 




















1 


V o 





1 


o J 



(2.15) 



where the permutation operator in coordinate representation is A = (0)(1, 2)(3). 
The operator U 2 has two eigenvalues (1 and -1), and the operator V 2 per- 
mutes between the vectors with the same eigenvalue (1 or -1). This is why 
the operator U 2 commutes with the operator V 2 . Permutation by the oper- 
ator A turns the operator U 2 into the operator AU 2 A\ which anticommutes 
with V 2 . The operators AU 2 A j; and V 2 form a complementary pair of oper- 
ators for sub-dimension Mi = 2, where their anticommutation is consistent 
with Eq. (11. 8p . The operators U 2 and AV 2 A^ form another complementary 
pair of operators for sub-dimension M2 = 2. The operator A permutes the 
eigenvalues of U 2 in such a way as to make the operator V 2 anticommute 
with AU 2 A\ 

A more interesting example is dimension M = 12, where both coprime 
and non-coprime factorizations are possible. For the case of M\ = 2 and 
M2 = 6, using the coordinate representation, the permutation operator is 



A=(0)(l,6,3,7,9,10,5,8,4,2)(ll). 
In the other case, where Mi = 3 and M2 = 4, the permutation operator is 
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A=(0)(l,4,5,9,3)(2,8,10,7,6)(ll). 

In both cases the construction of the operators from Eq. (12. 8p leads to the 
factorized pairs of operators. The generality of the new factorization enables 
to perform it for every factorized numbers Mi and M 2 . In the case where M\ 
or M 2 are composite numbers, another factorization can be performed until 
we reach prime numbers in factorization. 

Note that Schwinger's solution for non-coprime factorization in ref.[3j 
gives the factorized pairs of operators (see also ref.|13j). which obey the 
commutation relations of Eqs. ( II. 8| II. 9p . However, while for the coprime 
factorization an explicit expression is given in ref.[3] connecting between the 
factorized pairs and the original operators U and V, no such expression is 
given for the non-coprime case. In our paper this explicit expression is given 
in Eq. (jZED - 

3 New kq-\ike bases 

Each set (a') and (b) of operators (Eq. (I2.8j) ) generates M commuting op- 
erators and can be used for the definition of a basis for the M - dimensional 
Hilbert space. The set (b) of operators has, as an eigenbasis, the \K, Q) 
basis. As a result of the unitary transformation of the (a) set, the (a') set 

defines the kq-\ike basis \k,q) as follows: 



(a') \k,q) = A\k,q) = ^ ^TS' c*~|« + qM 2 ), 

(3-1) 

* = 5/, / = 0,...,M 2 -1, 
[a) \ q = 0,...,M 1 -l. 

As a result of the fact that the two sets of operators (Vi,Ui) and (V^, U2) 
describe Mi^and M 2 subspaces in the entire M - dimensional Hilbert space, 
the bases \k,q) and \K, Q) are mutually unbiased. Let us check the overlap 
between \k,q) and \K, Q) states (a = Mi, b = M 2 ): 

A/2-1 A/1-1 

(k,q\K,Q) = -j=^= Yl E e- lksa e lKtb (s + qM 2 \Q + tb). (3.2) 
Inserting (a = Mi, b = M 2 ) we have: 
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= Vm SSo" 1 ESo" 1 e~^e iKm ^s + qM 2 \Q + tM 2 ) = 
= EjlV 1 EST 1 e- ife ^ e ™^( S - Q)5 M2 (? - t) = (3.3) 

_ 1 c ~ikQM-i c iKgM?. 
\[M 

Here £ Ml (:r — Xo) means that the argument of the delta function is taken 
modulo Mi, 5 Mi (0) = 1 and elsewhere is zero. Therefore, these bases are 
mutually unbiased: | (k, q\K, Q) | 2 = -p. We call the basis \k,q) a kq-\ike basis 
because it is not an eigenfunction of the same operators as the \k,q) basis, 
but of the operators related to them by the permutation^transformation. In 
addition, it has different periodicity properties. As \k,q) is defined, it has 
the completely periodic property: 

(a') \k + ^-,Q) = \k,7+M 1 ) = \k~Tq). (3.4) 

To show explicitly the difference and similarity between the bases \k,q) 
and \k,q) we consider an example of dimension M = 6 with Mi = 3 and 
M 2 = 2. Using the \x) representation we list in three columns all basis 
members of \k,q) on the left hand side, all \k,q) basis vectors in the middle 
and \K, Q) on the right hand side: 



1(0,0)) = ^(|0) + |2) + |4)), 
|(0,1)) = ^(|1) + |3) + |5)), 
|(1,0)) = ^r|0) + e 2 f i |2) + e 4 ?|4)V 

|(l,l)> = ^Ul> + e^|3> + e¥|5>V 
l(2,0)) = ^(|0) + e ^|2)+ e a F|4)V 
l(2,l)) = ^(|l)+ e ^|3)+ e 2 i 1 |5)). 

(3.5) 

Hence, the \k, q) and the \k, q) bases are neither equal nor orthogonal to one 
another (they are eigenfunctions of different sets of operators) . Nevertheless, 
both these bases are mutually unbiased to the \K,Q) basis in this coprime 
case. 



' 1(0,0)) = ^(|0) + |3)) 
|(0,1)) = ^(|1) + |4)) 

1(0,2)) = ^(|2) + |5» 
' 1(1,0)) = ^(|0)-|3)) 
|(1,1)) = ^(|1)-|4)) 
|(1,2)) = -±=(|2)-|5)) 



' 1(0,0)) = 


7f (|0) + |1)), 


1(0,1)) = 


7f(l 2 ) + l 3 ))' 


1(0,2)) = 


^(|4) + |5)), < 


1(1,0)) = 


7f(|o>-|i», 


1(1,1)) = 


^ (|2> - |3>) , 


k 1(1,2)) = 


T7l(l 4 ) - 15)), 
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4 Engineering of the energy spectrum using the 
permutation operator A in Harper-like Hamil- 
tonians 

The new factorized pairs of operators (T'(a), t(^ l )) and (T(b),r , (— )) (Eq. 
(12. 8p ) describe M 2 and M\ dimensional subspaces, respectively [3j. Therefore, 
replacing the Harper-like Hamiltonian H[T(b), r(^)] of ref.[4j by H[T'(b),T(^)} 
we should obtain Mi energy levels, each of which is degenerate M 2 times, 
without restriction for M 2 and Mi to be coprime. 

To show the advantage of the new factorization, we compare it with the 
energy spectra design method of ref. [1]. As a first example, let us con- 
sider the dimension M = 6. We choose the simple Harper-like Hamiltonian 
proposed in ref. [I]: 

H = H(T(b), r(f )) = V x cos{\p) + V 2 cos(^x), (4.1) 

where V\ and V 2 are constants. We solve this Hamiltonian using the kq- 
representation: 

|V>) = IM(MV»> = 52C ktq \k,q). (4.2) 

k,q k,q 

The resulting eigenvalue equation for our Hamiltonian is: 

[V lC os{\p) + V 2 cos{^x)} E fe , 9 C k>q \k, q) = e £ M C k>q \k, q). (4.3) 
After applying the operators we have: 

E fc q CkM{\k, q-b) + \k,q + b)) + V 2 cos{^q)\k, q}} = S £ fe q C k>q \k, q). 

(4.4) 

The above eigenvalue equation can be solved for each value of k = jjf 
independently. So for our particular choice of dimension M = 6 with Mi = 
a = 2 and M 2 = b = 3, performing the summation over q values with the use 
of the quasi-periodicity property of the \k,q) states, we obtain the following 
equation (for some particular k value): 

C fc , [f ( e4fci l fc > !> + e ~ 2ki \ k > !» + ^cos(f • 0)|fc, 0)] + 
C k ^(e 2ki \k, 0) + e- 4fci |A:, 0)) + V 2 cos{\ -\)\k, 1)] = (4.5) 
= e[C kiO \k,0)+C k) i\k,l)}. 
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Using the orthogonality of the \k,q) states the above equation is equivalent 
to the solution of the following Mi = 2 coupled equations: 

[v ie ^f -v 2 )[c k , 1 )= £ [c k , 1 )- ^ 

The energy spectrum e is: 

e li2 = ± y/V? + Vi, (4.7) 

which is / independent and therefore each energy level is 3-fold degenerate 
(note that for current example / = {0,1,2} and k = ^-/). The relation 
between the coefficients is: 

y ie -¥ 4 / u ie ¥ 4 / 

C k ,o = 7J— C k,i or equally C M = — — C fcj o- (4.8) 

S — V 2 £ + V 2 

On the other hand, if instead of coprime factorized M = 6 we choose 
M = 4 and substitute Mi = a = 2 and M 2 = 6 = 2 into equation (I4.4p . we 
get Eq. (14. 9p with non-degenerate energy spectrum: 



E GUyd*. ff- 2 > + \ k > 9 + 2 )) + q)\k, q)]=sJ2 C k>g \k, q). (4.9) 

k,q k,q 

Using the quasi-periodicity properties of \k, q) states we have: 



£ C k , q [V lC os(2k)\k, q) + V 2 cos(—q)\k, q)) = e ^ C k>q \k, q), (4.10) 

fe,g k,q 

and the energy spectrum is: 

£1,2,3,4 = ±V X ± V 2 . (4.11) 

This result is expected, because of the absence of factorization into sub- 
dimensions Mi = 2 and M 2 = 2 using the operators of Eq. (II. lip . 

Let us now follow the same procedure with the new operators of Eq. (12. 8p . 
Accordingly, the Harper-like Hamiltonian of Eq. (14. ip changes to: 

H = F(T(6),r'(f )) = V x cos{\p) + V 2 Acos(^x)Al (4.12) 

To compare the two schemes we solve the above Hamiltonian using the kq- 
representation: 
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fc,<j fc,q 

The eigenvalue equation for our Hamiltonian is: 

[ Vl cos{\p) + y 2 Aco S (^x)At] E fc , ff ^IM> = eE k , g C k , q \Kq)- (4.14) 
After applying the operators we have: 

(4.15) 

where we have used the two relations: 

T{b)\k^q) = %q^\) and r'(— )|M) = e^M)- 

a 

As before, the eigenvalue equation (14. 15ft can be solved for each value of k 
independently, and using the complete periodicity property of the \k,q) we 
have (with Mi = a = 2 and M 2 = 6 = 3): 

<7 fc ,o[£(|M) + IM)) + • 0)|M>]+ 

c M [f(M> + |M» + ^Wf • 1)|M)] = (4-16) 

= e[C fc ,o|M> + C r MlM)]- 
In matrix form the above equation reads: 

(S-'U (&)-■(&)■ 

Hence, we get the same spectrum of energies as before, with each level being 
3- fold degenerate: 

e 1)2 = ±y/V? + VI (4.18) 

and a new relation between the coefficients: 

Vx - Vx ~ 

C k ,o = 77CV1 or equally C k) i = — —rrC kfi . (4.19) 

E — V2 6 + V2 

In comparison with the previous solution Eq. ()4.8|1 . now we do not have a 
k-dependent phase in the relation between the coefficients C k ,o and C k) \. 
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In the case of M — 4, solving equation (14.15)) with Mi = a = 2 and 
M 2 = b = 2, we have: 



^^.JydM^l) + 1» + V 2 cos(^q)\k,q)] =eJ2C k , q \k,q). 

k,q k,q 

(4.20) 

As a result of the k independence of the equation above it is equivalent to 
the eigenvalue equation considered for dimension M — 6. Therefore (as one 
can easily check) we have to solve the matrix equation: 




(4.21) 



and consequently the corresponding energy levels, with each level being 2-fold 
degenerate, are: 

e 1)2 = ±y/V 1 2 + V 2 2 . (4.22) 

Therefore, in spite of the non-degenerate spectrum of the Hamiltonian H(T(b), 
for non-coprime Mi and M 2 , for the Hamiltonian H(T(b), r'(— )) the energy 
levels preserve their degeneracies. 



5 Summary and discussion 

The main result of our work is a generalization of the Schwinger unitary 
operator factorization to non-coprime factorizations. That is, for a composite 
dimension M = MiM 2 , we factorize the U and V operators from Eq. (11. 3p 
into two pairs of operators (Ui,Vi) and (U 2 ,V 2 ) Eq. (12. 8ft . Each of the 
pairs generates a complete orthogonal operator basis for the sub-dimensions 
Mi and M 2 , and operators from different bases commute. The factorization 
enables us to consider any single physical system with dimension M = MiM 2 
as a pair of physical systems in Mi and M 2 - factorized degrees of freedom, 
where Mi and M 2 are not restricted to be coprime. Considering factorized 
operators may simplify various M - dimensional phase space problems in the 
same way as the Cooley-Tukey FFT simplifies the application of the DFT. 
Moreover, the new factorization deepens our physical intuition. In particular, 
we applied the new factorization to a Harper-like Hamiltonian model, and 
developed an algorithm for energy spectrum design in this model. Using 
the algorithm, we can construct a Hamiltonian, which is a function of the 
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operators (U\, V\). Therefore, it is designed to obtain M± energy levels (with 
a spectrum determined by Hamiltonian's details), each level being M 2 -fold 
degenerate. The algorithm of the energy spectrum design can be of interest, 
for example in solid state physics for electrons in a strong magnetic field [14j. 

The application of the permutation operator A (which is the key to the 
solution for the non-coprime cases) to the kq bases problem generates the kq- 
like basis which is a MUB to the original \K, Q) basis. This kq-\ike basis has 
a different periodicity property than the original kq bases: it is completely 
periodic in the coordinate and momentum variables simultaneously. 
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